Frustration induced Raman scattering in CuGeOa 



< 

m 



> 
O 
00 

O 



o 



V. N. Muthukumar^, Claudius Gros^, Wolfgang Wenzel"'^, Roser Valenti^, 
P. Lemmens^, B. Eisener^, G. Giintherodt^, M. Weiden^, C. GeibeP, F. Steglich^ 
^ Institut fiir Physik, Universitdt Dortmund, 44^'^^ Dortmund, Germany 
^ 2. Physikalisches Institut, RWTH Aachen, 52056 Aachen, Germany 
^ FB Technische Physik, TH-Darmstadt, Hochschulstr. 8, 64289 Darmstadt 

(February 1, 2008) 

We present experimental data for the Raman intensity in the spin-Peierls compound CuGeOs 
and theoretical calculations from a one-dimensional frustrated spin model. The theory is based on 
(a) exact diagonalization and (b) a recently developed solitonic mean field theory. We find good 
agreement between the ID-theory in the homogeneous phase and evidence for a novel dimerization 
of the Raman operator in the spin-Peierls state. Finally we present evidence for a coupling between 
the interchain exchange, the spin-Peierls order parameter and the magnetic excitations along the 
chains. 
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Low-dimensional spin systems exhibit many unusual 
properties resulting from quantum and dimensionality 
effects. An example is the continuum of spin- wave ex- 
citations in quantum onc-dimensional (ID) spin systems 
which has been predicted for a long time and has re- 
cently been confirmed by neutron scattering experiments 
§ on KCuFg. 

Quantum ID spin systems also show a variety of in- 
stabilities. Of particular interest is the spin-Peierls (SP) 
phase due to residual magnetoelastic couplings , which 
leads to the opening of a gap in the spin excitation spec- 
trum. The discovery [Q of the spin-Peierls transition be- 
low Tgp 14 K in an inorganic compound, CuGeOa, has 
attracted widespread attention. This compound consists 
of chains of spin- 1/2 Cu^+ ions coupled by antiferromag- 
netic superexchange via the oxygen orbitals ||J^ . The Cu 
ions lie along the crystallographic c-axis and the exchange 
along the chains can be modeled by the ID Hamiltonian 



H = J^[(l + ,5(-l)^)S, • S,+i + aS, • S,+2] 



(1) 



where 5 is the dimerization parameter that vanishes 
above Tsp . The special geometry [^,|| of the super- 
exchange path in CuGeOs leads to a small value of the 
exchange integral J w 150K and a substantial n.n.n. 
frustration term ~ a which competes with the n.n. anti- 
ferromagnetic exchange. The interchain couplings have 
been estimated to be small, Jb ~ O.IJ and Ja ~ — O.OIJ 
for the interchain exchange constants along a- and b- di- 
rections, respectively 

The phase diagram of H in Eq. (|l]) has been calculated 
using the density-matrix renormalization-group method 
1^. For 5 — and a < ac ~ 0.2411, the ground state is 
gapless and renormalizes to the Heisenberg fixed point. 
For a = 0.5 and i5 = 0, the ground state is given by a 
valence-bond state and a gap of order J/2 induced by 
frustration is present. While the evaluation of the dy- 
namical properties of ([^) is a challenge to theory, the 
rich phase diagram can be explored by a variety of in- 



teresting experiments. In this context, the substantial 
value of the n.n.n. exchange integral in CuGeOa allows 
the experimental investigation of the effects of compet- 
ing interactions in a low dimensional magnet, both in the 
uniform and in the spin-Peierls state. 

An experimental method particularly suited for the 
study of magnetic excitations in an antiferromagnet is 
two magnon Raman scattering. For CuGeOs, the Ra- 
man operator in Aig symmetry [Qis proportional to 



-7(-l)^)S.-S 



i+l 



(2) 



In the homogeneous state ((5 = 7 = 0) the interaction 
Hamiltonian commutes with the Heisenberg Hamiltonian 
for the case a = and there would be no Raman scatter- 
ing. However when a ^ 0, the model (|l|) leads to mag- 
netic Raman scattering due to the presence of compet- 
ing interactions which can be observed experimentally. 
Previous fits to experiments jll) based on ([|) assumed 
a = 0. Next, we note the presence of the factor 7 in 
Eq. (H) which appears because the exchange integral is 
sensitive to the inter-ionic distance. 

Both features mentioned above are taken into consid- 
eration in this Letter where we present experimental data 
as well as the first theoretical results for the Raman in- 
tensity in one-dimensional spin systems obtained from 
(a) exact diagonalization studies of chains with Ng <28 
sites and (b) a newly developed solitonic mean field the- 
ory. Our results for T > Tsp are obtained with (5 = 
in Eq. (|^), where the mean field theory provides an an- 
alytical framework to obtain expressions for the Raman 
intensity at finite temperatures. Below Tsp, we present 
first experimental evidence for a qualitative change in 
the Raman operator at the SP transition through the 
appearance of the factor 7 in Eq. (|^). Here, we do not 
attempt to calculate 7 microscopically, but by compar- 
ing with experiment. We find that 7 « 0.12 and this 
value of 7 affects the shape of the Raman spectrum in 
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the SP phase dramatically. Indeed, we find that the ex- 
perimental results for T < Tsp cannot be explained by 
Eq. (|l|) assuming 7 to be zero. A detailed comparison 
of experiment and theory, both below and above the SP 
transition, is presented. 

We first give a brief account of the theoretical and 
experimental methods before discussing the experimental 
data. The Raman intensity at zero temperature is given 
by 

/flH - --/m(0|ilfl— . l—-—Hn\0), (3) 

TT uj + le — (H — La) 

where Eq is the ground-state energy, H the Hamiltonian 
given by Eq. (|l|) and e —>■ 0+. For a numerical evalua- 
tion of Eq. (|), we scale the Hamiltonian H = cX + d 
such that the eigenvalues of the rescaled Hamiltonian 
X are in the interval [—1,1]. We define a rescaled en- 
ergy and frequency by Eq = cxq + d and w = cx + d 
and expand Ir{x) in terms of Tschebycheff polynomi- 
als, Ti{x): Ir{x) = [l - x^)-^/^ Yl=oaiTi{x + xq) with 
a; = (2 - 5i,q)/^{Q\HrTi{X)Hr\Q) and Np ^ 00. The 
quantities a/ are evaluated recursively. This procedure, 
the kernel polynomial approximation, is an established 
1^ and numerically stable method for the evaluation of 
the density of states For large values of iVp, the re- 
sulting spectral weight consists of a series of very sharp 
peaks that become delta- functions in the limit Np — > 00. 
Generally we find Np — 100 convenient for comparison 
with experiment. 

We also calculate the Raman intensity above Tsp us- 
ing a solitonic mean field theory. This method, orig- 
inally suggested for the n.n. Heisenberg chain |Q, is 
based on a transformation of the Heisenberg Hamiltonian 
to a Hamiltonian describing the dynamics of antiferro- 
magnetic domain walls (solitons or spinon excitations) 
and a subsequent Jordan- Wigner transformation. The 
method reproduces the exact solutions at both the Ising 
and XY-limits. For the Heisenberg model it leads to cor- 
rect asymptotic behavior of dynamical correlation func- 
tions. The solitonic mean field theory leads to a ground 
state described by the Hamiltonian Hs = J2k -^kcj^Ck, 
where the c^'s are quasiparticle (spinon) operators that 
are linear combinations of the soliton operators. The 
mean field dispersion relation Ek — {1 + ^)J\ cosfcaj com- 
pares very well with the des Cloizeaux-Pearson spectrum 
^Jjcosfcaj obtained from Bethe ansatz |l^. We gener- 
alize this method to the case of the n.n.n. chain. 

The absence of Raman scattering in the n.n. chain 
is understood as Hs conserving the number of spinous. 
This picture changes when a ^ 0. We find that the inclu- 
sion of the n.n.n. term generates the following processes: 
(i) two-spinon scattering terms that lead to a renormal- 
ization of the spinon velocity and (ii) four-spinon cre- 
ation terms that generate two-magnon Raman scattering 
in the n.n.n. chain. A perturbative treatment of the 



four-spinon terms allows us to obtain expressions for the 
Raman intensity Ir{lo) in the homogeneous phase at all 
temperatures. 

The experiments were performed using the excitation 
line A = 514.5-nm of an Ar-laser with a laser power of 
2.7mW. We ensured that the incident radiation does not 
increase the temperature of the sample by more than 1.5 
K. We used a DILOR-XY spectrometer and a nitrogen 
cooled CCD (back illuminated) as a detector in a quasi 
backscattering geometry with the polarization of incident 
and scattered light parallel to the c-axis and the Cu-0 
chains, respectively. Details of the experiment will be 
published elsewhere ||l6|] . 

In Fig. |l|, we present the data for the two-magnon Ra- 
man continuum in the homogeneous state at T = 20 K. 
Phonon lines ||ll| at 184 cm~^ and at 330 cm^^ are sub- 
tracted from the experimental data (squares). The shoul- 
der observed at ~ 390cm^^ is presently not understood. 
We interpret the two-magnon continuum as scattering in- 
tensity caused by the creation of four spinon excitations. 
The Raman intensity from such a scattering process cal- 
culated from our solitonic mean field theory is shown in 
Fig. |l] (dashed-dotted line), where we use J = 150 K ~ 
104 cm~^ and a = 0.24 The maximum theoretical 
value of the Raman intensity is normalized to the exper- 
imental value from which a Rayleigh tail and a uniform 
background of 50 counts are subtracted. The maximum 
of the experimental data is situated at slightly larger en- 
ergies, which could be a consequence of a slight misfit of 
the parameters used. For instance, we find that J = 160 
K and a = 0.2 (dotted line) improves the agreement be- 
tween theory and experiment. The data presented in Fig. 
|l| provide substantial evidence that a ID Hamiltonian of 
type (P can indeed account for the observed Raman con- 
tinuum owing to the presence of the n.n.n. frustration 
term. This continuum should also be seen in neutron 
scatteringexperiments above Tsp. 

In Fig. y we present the numerical results for the dimer- 
ized state, S = 0.03 and a finite-size analysis (in the 
inset). The numerical result for 7 = (upper curve, 
dashed line) is surprisingly flat and can be approximated 
(upper curve, solid line) by the expression 

Ir{j = 0, w) = Ae{Lj - 2A) (1 - tanh[2(w - ujq)]) , (4) 

with the values of the parameters 2 A « 30cm~^, ujq « 
312cm~^ being determined by a fit to the numerical data 
(A is a normalization constant). Below Tsp, a gap opens 
up and the DOS has a singularity at the lower edge 
|ll|,0- Therefore, one expects a peak arising from this 
singularity. However, this singularity is removed by ma- 
trix element effects arising from the Raman operator [ p^ . 

The 7-dependence of the Raman intensity is strong, as 
can be seen in Fig. |^, where we plot the data for 7 = 0.12 
(lower curves) In order to understand this large matrix- 
element effect we have examined in detail (for Np — 1000) 
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the relative weight ^(7, Ei) of the individual poles (at en- 
ergies Ei) contributing to the Raman intensity for sys- 
tems of size Ng = 20, 24 and 28. We have found that the 
basis of the observed matrix-element effect lies in the re- 
markable fact that for each pole there is a certain jo{Ei) 
at which the intensity of the pole actually vanishes, i.e. 



7 - lojEj) 



(5) 



where Ij is a normalization constant defined by 
/p " dujp{j,uj) — 1 (we choose ujc = 6J). From the nu- 
merical data wc found that 70 (-Ei) can be approximated 
by 7o(£'i) ~ (5 + nEi, where k w 4/3000cm. Combining 
(^ and (^, we obtain an analytic approximation for the 
Raman intensity in the dimerized state, 



(6) 



This formula reproduces the numerical results well (see 
Fig. I). 

In Fig. we present data for the spin-Peierls phase at 
T = 5K. Let us first discuss the line at 30 cm~^. We find 
that the analytic curve for 7 = 0.12 reproduces this peak 
well (solid line in Fig. Note that here 7 is the only free 
parameter in the theory. Choosing 7 = instead would 
result in a complete disagreement with experiment (com- 
pare Fig. ^). The line at 30 cm^^ is known to be of ID 
magnetic origin (the value 30 cm~^ is indeed twice the 
one-magnon gap obtained from neutron scattering [^). 
This suggests that this peak arises from the spin-Peierls 
Hamiltonian (|l|) in combination with matrix-element ef- 
fects of the Raman-operator (^ . Therefore we conclude 
that Fig. ^ provides strong evidence for the dimerization 
of Hn below Tsp leading to substantial matrix-element 
effects. 

Let us now consider the two lines at 226 cm^^ and 
at 104 cm~^ observed below Tsp (see Fig. ||). The as- 
signment of the 104 cm~^ line is still controversial and 
will not be discussed here. The 226 cm^^ line has been 
assigned previously to be of magnetic origin |Tl| ] . A clas- 
sical (non-interacting) spin-wave calculation, using the 
measured magnon dispersion in c-direction 1^ produces 
a peak at 226 cm~^ jl^. However, it is well known that 
in ideal one-dimensional systems, magnons do not behave 
classically but form a continuum of excitations Q . Con- 
sequently one would not expect a sharp peak, like the 
one observed at 226 cm~^ to occur in a one-dimensional 
system. This expectation is borne out by our numerical 
data (solid line in Fig. |3|), which does not show the 226 
cm~^ peak. It is clear that to explain the presence of 
this peak, one has to go beyond a one-dimensional spin 
model. 

The magnitude of the interchain coupling Jh in 
CuGeOa is comparable to the spin-Peierls temperature 
(14 K). One might therefore expect the chains to become 
correlated for temperatures around Tsp « Jb, resulting 



in the appearance of a well defined magnon branch, which 
leads to the observed peak at 226 cm~^. In Fig. ^, we 
present the temperature dependence of the intensity of 
the 226 cm~^ line (for comparison we include in Fig. |^, 
two other Peierls- active Raman lines). One clearly sees 
that this line becomes active only below Tsp- Thus we 
conclude that the interchain coupling becomes relevant 
only below Tsp- This conclusion is in agreement with 
recently presented neutron scattering data |jl^, which 
shows a pronounced change in the spectrum below Tsp. 

In conclusion, we find good agreement between the 
experimental magnetic Raman spectrum in the homo- 
geneous phase and the one-dimensional frustrated spin 
model for CuGeOs. The importance of the competing 
interactions for the occurrence of magnetic Raman inten- 
sity in the homogeneous phase has been pointed out and 
is consistent with the absence of inelastic Raman scat- 
tering in ID spin compounds without frustration such 
as KCuFs 11^]. For the spin-Peierls state we find that 
observation of a sharp line at the spin-triplet excitation 
energy of 30 cm^^ indicates a dimerization of the Raman 
operator. Finally, the observed Peierls-active line at 226 
cm~^ indicates the appearance of a well defined magnon 
branch below Tsp- 
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FIG. 1. Experimental (squares) and theoretical results 
(lines) for the Raman intensity in the homogeneous phase ol 
CuGeOa at T = 20 K. Shown are the results from the solitonic 
mean field theory at T = 20 K (5 = = 7) and for two sets 
of parameters (a) J = 104 cm~^, a — 0.24 M (dashed-dotted 
line), and (b) J = 119 cm"\ a = 0.2 (dotted line). 



FIG. 2. The numerical results for Ns = 28 (dashed lines), 
5 = 0.03, a = 0.24, Np = 100 and 7 = (upper curve) and 
7 = 0.12 (lower curve). The numerical results are compared 
with the analytic formula ^ (solid lines). Inset: A com- 
parison of numerical results for Ns — 28 (dashed line) and 
Ns = 24 (dotted line) and 7 = 0. 

FIG. 3. Experimental data (squares) and the analytic ap- 
proximation (Eq. ^ solid line) for the Raman intensity of 
CuGeOa in the spin-Peierls phase at T = 5 K. The phonon 
lines have been subtracted from the experimental data. 

FIG. 4. Temperature dependence of the relative intensi- 
ties of three Peierls-active Raman lines. The 30 cm~^ line 
(squares), the 226 cm"'^ line (triangles) and the 369 cm"'^ 
line (circles, a folded phonon line) The intensities are normal- 
ized to their respective values at T = 6 K. The intensities at 
T = 200 K are subtracted for reference. 
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